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We have analysed the first 15 months of Fermi/LAT data of the radio loud quasar 3C 273. Intense gamma- 
ray activity has been detected, showing an average flux of F{> 100 MeV) = 1.4 X 10"^ ph cm s , with a 
peak at F(> 100 MeV) = 5.6 X 10"'' phcm'^ g-l detected during a flare in September 2009. Together with 
the brightening of the source, a possible hardening of the gamma-ray spectrum is observed, pointing to a 
shift of the inverse Compton peak toward higher energies than the 1—10 MeV range in which 3C 273 inverse 
Compton emission is typically observed to peak. During the 15 months of observations the photon index is 
measured to vary between F = 2.4 and F = 3.3, with an average value of (F) = 2.78 ± 0.03. When compared 
to the observations at other wavelengths, the gamma-rays show the largest flux variations and we discuss the 
possibility that two different components are responsible for the inverse Compton hump emission below and 
above the MeV peak. 



1. Introduction 



3C 273 is a bright, radio loud quasar that shares 
most of the characteristic properties of blazars, Uke 
strong radio emission, a jet with apparent superlu- 
minal motion, large flux variations, and, occasionally, 
polarisation of the optical emission (see Courvoisier et 
al. 1998 for a review). On the other hand, this AGN 
presents also characteristics more common in Seyfert 
galaxies, like a prominent blue bump Q and the signa- 
ture of non-thermal Comptonisation processes in the 
X-ray spectrum 3]. 

3C 273 was the first extragalactic source to be 
detected at gamma-rays by COS-B observations in 
1976 |3| and later on was detected several times 
by the EGRET instrument on board the CGRO 
satellite [5]. Thanks to its early discovery as 
a multiwavelength emitter and to its brightness, 
3C 273 is one of the best observed AGN at all 
energies. The multiwavelength on-line database of 
3C 273 {http://isdc.unige.ch/3c273/) contains about 
40 years of data from radio to gamma-rays [6|. This 
set of data constitutes a very precious starting point 
for studying the multi-band variability of this ob- 
ject. In particular, with the launch of Fermi in June 
2008 an intense gamma-ray monitoring of the source 
has started and triggered observations along the 
whole electromagnetic spectrum, coordinated by sev- 
eral groups. This extensive and regular look at the 
source, together with the comparison with the histor- 
ical long-term behaviour at all wavelengths, will allow 
us to investigate the properties and origin of the high- 
energy emission in greater detail than it had ever been 
done before. 



2. From radio to X-rays 

2.1. The 3C 273 multiwavelength 
database 

The 3C 273 database, first created in 1999 [1], 
has been recently updated with observations up to 
2007 [G] in order to provide the community with one 
of the most complete multiwavelength databases for 
an AGN. At the time being the database contains 
70 light curves covering 17 decades in energy and 
spanning up to 40 years of observations (Fig. [J). 
In addition, for the X-ray band, information about 
the spectral parameters used to extract the fluxes 
is also available. The bulk of the high-energy data 
has been collected in the last 10-15 years, thanks 
to the intense X-ray monitoring provided by Bep- 
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Figure 1: Time coverage of all 3C 273 observations of the 
database as a function of the frequency, color-coded to 
represent the vFi, intensity of the measured fluxes. 
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Figure 2: Spectrum of the fractional variability 
amplitude, computed using the light curves of the 3C 273 
database. Above 100 MeV, the circles represent the 
variability observed with EGRET data and the open 
squares are the results obtained by the LAT. 



poSAX, CGi?0/BATSE, RXTE, XMM-Newton, IN- 
TEGRAL and Swift. In the gamma-ray domain, 
the data collected by CGi?0/OSSE and EGRET 
have already been included in the database, while 
we plan to include also results from AGILE /GKID 
and Fermi /hAT observations, within the upcoming 
months. 

2.2. The variability study 

As a first application of the on-line database, we 
studied the variability properties of 3C 273 at differ- 
ent wavelengths and, correlating the behaviour in the 
different energy bands, we tried to identify the emis- 
sion processes and understand their inter-connections 

The amplitude and the time scales of the variations 
of 3C 273 depend strongly on the frequency and show 
trends that are characteristic of the underlying emis- 
sion processes (Fig. [2|). For example, the increasing 
amplitude of the variations and the decreasing time 
scales when going from long radio to short millimeter 
wavelengths is the typical signature of the superposi- 
tion of synchrotron flares from shock waves propagat- 
ing along the jet. 

The different variability properties in the X-rays be- 
low and above 20 keV seem to support the possibility 
of two separated emitting components, possibly non 
thermal Comptonisation (as in Seyfert galaxies) and 
inverse Compton from the jet (as in blazars). As an al- 
ternative, a single component characterised by at least 
two independently varying parameters could also ex- 
plain the X-ray timing properties. 

In any case, the dominant hard X-ray (> 20keV) 
emission is most probably not due to electrons accel- 



Figure 3: Correlation between the hard X-ray (20-70 
keV) and the 1.3 mm light curves. No correlation is 
detectable at short time lags. 

erated by the shock waves in the jet as their variability 
does not correlate with the flaring millimeter emission 
at short time lags (Fig. [3]). Instead, the hard X-rays 
seem to be associated to long-timescale variations in 
the optical domain. This optical component is consis- 
tent with being optically thin synchrotron radiation 
from the base of the jet and the hard X-rays would 
be produced through inverse Compton processes (syn- 
chrotron self-Compton and/or external Compton) by 
the same electron population. 

In the gamma-ray domain, the limited data avail- 
able before AGILE and Fermi observations did not 
allow to draw a firm conclusion about the variability 
in the MeV-GeV energy range. The increase of the 
variability amplitude observed in the X-rays seems to 
continue up to the gamma-rays, possibly even above 
100 MeV where the observed decrease is likely due to 
the large uncertainties of the EGRET fluxes (circles 
in Fig.[2|). 



3. The gamma-ray emission 
3.1. Before Fermi 

During CGi?0/EGRET observations, 3C 273 was 
either found below the EGRET sensitivity level, or de- 
tected at a relatively low gamma-ray flux, with an av- 
erage fluxof F(> 100 MeV) = 15 x 10"** phcm^^ g-i, 
between 1991 and 1995, and a maximum detected flux 
of F(> lOOMeV) ~ 80x 10-**phcm-2s-i in 1997 ^. 
The EGRET average spectrum has a power law shape 
with photon index (F) = 2.4 ± 0.1, but variations of 
the spectral slope have been observed, with the spec- 
trum hardening while the source brightened ^, il] ■ A 
flux variability by a factor of 4 within one week is 
the most extreme behaviour observed for 3C 273 at 
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Figure 4: Fermi /LAT count map of the 15° region 
around 3C 273, produced with data collected during 15 
months. Two additional sources are clearly detected: the 
blazar 3C 279 and the unidentified source at 
R.A. = 189.9, DEC = 4.7. 

gamma-rays by EGRET 0. Multiwavelength cam- 
paigns have been carried out simultaneously to the 
the gamma-ray observations in 1993-1995 and showed 
that the emission detected by EGRET can be equally 
well explained with synchrotron self-Compton (SSC), 
external Compton (EC) or proton-initiated cascade 
(PIC) models 

During the first year of observations, A GILE / GRID 
has detected 3C 273 with an average flux of F(> 
100 MeV) = 24 X IQ-Sphcm-^ s"\ consistent with 
EGRET measurements [ifl]- During the multiwave- 
length campaign of December-January 2008, the 
source was detected only during the second week of 
observations at a flux of F(> 100 MeV) = 34 x 
10^^ phcm~^ [ll|. The hard X-ray to gamma-ray 
emission is modeled as external Compton radiation on 
thermal disc photons, and the variability of this high- 
energy component is consistent with an acceleration 
phase of the electron distribution that moves the in- 
verse Compton peak to higher energies. 

3.2. Fermi/LfiJ data analysis 

The Fermi satellite, launched in June 2008, car- 
ries two instruments able to observe the gamma-ray 
sky: the Large Area Telescope (LAT), a pair con- 
version gamma-ray telescope sensitive to gamma-rays 
from 30 MeV to 200 GeV [ll; and the GLAST Burst 
Monitor (GBM), observing the sky in the 8 keV to 
30 MeV energy range |13|. The main instrument 
LAT is 10 times more sensitive than its predecessors 
CGi?0/EGRET and AGILE/GRID and, thanks to 



its large field of view of ~ 2.4 sr, it can scan the full 
sky every 3 hours. Already after the first 3 months of 
observations, the LAT had detected 205 objects with 
high-confidence, of which 106 are AGN 0,[il, includ- 
ing 3C 273. 

We analysed the Fermi /LAT data collected dur- 
ing the first 15 months of observations, between Au- 
gust 10, 2008 and October 19, 2009. We first down- 
loaded the 0.1-300 GeV light curve of 3C 273 pro- 
vided through the Fermi Monitored Source List Light 
Curves page^, selected only the detections on daily 
basis and rebinned the light curve to 3-day bins. We 
then performed the analysis on each of these 3-day 
time intervals selecting the event data with the gts- 
elect tool distributed as part of the Fermi-LAT Sci- 
enceTools (version 9.15.2). For each bin (correspond- 
ing in average to an integration time slightly smaller 
than 200 ks), the count map and exposure map were 
generated. 

For source detection and flux estimate, we used 
the gtlike tool, based on a maximum likelihood al- 
gorithm. For a region of interest (ROI) of 15° around 
3C 273, a model containing the isotropic and Galactic 
diffuse background models (as given by default) and 
all the sources in the field was fitted to the data. For 
the field and data analysed, beside 3C 273 two addi- 
tional sources were considered: the blazar 3C 279 and 
the unidentified source at R.A. = 189.9, DEC = 4.7 
(Fig.d]). All the sources have been modeled with sim- 
ple power laws in the energy range 0.1-200 GeV, in or- 
der to obtain their spectral parameters and test statis- 
tic, TS. We note that for each source the square root of 
TS provides an approximate estimate of the detection 
significance. 

Following 0], we added 3% systematics to the un- 
certainties on the estimated fluxes. 

3.3. Fermi /LAT preliminary results 

Figure [5] shows the results obtained during the first 
15 months of observations. In the selected time in- 
tervals analysed, 3C 273 was always detected with 
VTS > 5. Since the launch of Fermi, 3C 273 has 
undergone several episodes of very intense gamma- 
ray activity. The average flux above 100 MeV, 
F{> 100 MeV) = 1.4 x 10"*^ phcm'^ s"^ = 5.2 x 
2^0-10 ergcm"^ s^^, is significantly larger than the his- 
torical measurements of EGRET, with a peak flux of 
F{> 100 MeV) = 5.6 X 10-6 phcm-^s-i detected dur- 
ing a large flare in September 2009. The amplitude of 
the variations, Fvar (quantified through the fractional 
variability amplitude; see Q for more detail), mea- 
sured by Fermi, is consistent with the EGRET one at 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/msLlc/ 
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Figure 5: LAT results after about one year of 
observations of 3C 273. From the top to bottom panels: 
square root of the test statistic (corresponding to the 
detection significance); integration time for each data 
bin; gamma-ray flux above 100 MeV; photon index from 
a simple power law fit. In the bottom panel, the line 
indicates the average photon index over 15 months and 
the shaded area indicates the dispersion from the mean. 
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Figure 6: Zoom of Figure [S] on the last 3 months of 
observations. During the large flare of September 2009, 
the spectrum of the source seems to be significantly 
harder then outside the flaring period. 



brighter" is clearly observed (Fig. [7|): when 3C 273 
was brighter than 2 x 10~^ phcm~^ s~^, the photon 
index was always smaller than 2.8. The same effect 
was observed also in EGRET data during two flares 
of 3C 273 @,i and seems to be a common behaviour 
of the gamma-ray emission of blazars (see for example 
3C 454.3 and AO 0235-^164; [H, HI]). As proposed 
by Pacciani et al. (2009) to explain the gamma-ray 
variability observed with AGILE, the inverse Comp- 
ton peak might be shifted toward higher frequencies 
following an accelaration episode of the electron pop- 
ulation. Since the inverse Compton peak of 3C 273 
usually falls in the 1-10 MeV range jl, [T^] (Fig. E]), 
a shift toward higher energies could indeed produce a 
harder spectrum in the LAT energy band. 



100 MeV, but it is significantly larger above 100 MeV, 
confirming the increasing trend of i^var already ob- 
served in the X-ray domain (Fig. [2]). The decrease of 
i^var observed with EGRET above 100 MeV was likely 
due to its limited sensitivity. In addition, the ampli- 
tude of the variations measured in the gamma-rays by 
Fermi is the highest across the whole electromagnetic 
spectrum, confirming the extreme variability nature 
of the gamma-ray emission. 

The bottom panel of Figure [5] shows the evolution 
of the gamma-ray photon index during the LAT ob- 
servations. The photon index varies between F — 
2.4 ± 0.2 and F — 3.3 ± 0.3, with an average value 
of (F) = 2.78 ± 0.03. In general, the variations of 
the spectral slope do not seem to be significant, when 
the uncertainties are taken into account, with the ex- 
ception of a decrease of the photon index when the 
source brightened during the last 3 months of obser- 
vations (Fig. [HI). In fact, when plotting the photon 
index versus flux, the source behaviour "harder when 



4. The inverse Compton peak 

Gamma-ray observations of 3C 273 are essential to 
constrain the overall shape of the inverse Compton 
branch and to better understand the connection be- 
tween the jet emission in different energy domains. 
The most recent radio to gamma-ray campaigns car- 
ried out on this object, together with the numerous 
ones performed in the past, provide a large set of data 
to model the spectral energy distribution (SED) dur- 
ing several epochs and with different integration times 
(the average SED in Fig. [H). 

Multiwavelength campaigns performed during 
EGRET observations could not discriminate between 
the most common models used to describe the high- 
energy emission of blazars [5j] . A different variability 
behaviour below and above ^ 3 MeV suggested 
that two different mechanisms are responsible for the 
emission in the MeV and in the EGRET domain, with 
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the latter one being most likely produced by inverse 
Compton scattering of accretion disc photons [s]. On 
the other hand, the different variability of the hard 
X-rays and gamma-rays detected by INTEGRAL 
and AGILE has been explained by a single external 
Compton component shifted toward higher energies 
following an acceleration of the electron population, 
whereas a synchrotron self-Compton contribution is 
evident at lower energies [ll| . 

Several multiwavelength monitoring programs are 
on-going to follow the evolution of 3C 273 and many 
other gamma-ray sources. The preliminary results 
presented during the 2009 Fermi Symposium for 
3C 273 indicate that there is a correlation between the 
gamma-ray flaring states and the brightening of the 
VLBI core, with gamma-ray flare following the ejec- 
tion of a superluminal knot by up to 70 days [l^ . This 
would indicate that the gamma-ray high states are re- 
lated to a disturbance crossing the millimeter-jet core. 
An increase of the radio emission at 15 GHz has been 
observed quasi-simultaneously to the large September 
2009 gamma-ray flare, and the optical flux started also 
to increase in July/ August 2009, but no data were 
available at the time of the gamma-ray flare [l^. In 
addition, a correlation between the small gamma-ray 
flare of April 2009 and the optical emission is clearly 
detected p'l- 

As the flaring optical emission of 3C 273 is believed 
to be synchrotron emission related to the millime- 
ter jet, a gamma-ray/optical correlation would sup- 
port the origin of the high-energy emission in the mil- 
limeter jet and could indicate that the same electron 
population is responsible for both emissions through 
synchrotron and inverse Compton processes. On the 
other hand, we found the hard X-ray emission to be re- 




Figure 7: LAT gamma-ray photon index versus 
gamma-ray flux of 3C 273. Blue (black) points represent 
detections with TS larger (smaller) than 100. The source 
behaviour "harder when brighter" is observed in the tail 
at fluxes above 2 x 10~® phcm~^ s~^ with photon indices 
always below the average value of F = 2.8. 
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Figure 8: Average spectral energy distribution of 3C 273 
over up to 40 years. The error bars represent the 
standard deviation from the mean values and the green 
area indicates the observed range of variations. Above 
30 MeV, the circles correspond to the average EGRET 
data, whereas the blue cross hatched area indicates the 
flux range of 3C 273 during Fermi /LAT observations. 



lated rather to another optical component, varying on 
longer time scales, and produced as synchrotron emis- 
sion further away from the millimeter jet 0]. There- 
fore, this would point to a different origin of the hard 
X-ray and gamma-ray emissions. 



5. Conclusions 

Thanks to its large field of view and its unprece- 
dented sensitivity above 100 MeV, the LAT instru- 
ment is providing us with long-term, well sampled 
light curves and spectral information for a large num- 
ber of AGN [il [IE [in. Among them, 3C 273 
is a bright quasar extensively observed at all wave- 
lengths and showing a particularly intense gamma- 
ray activity since the launch of Fermi. This source 
was detected during the first 15 months of obser- 
vations with an average fiux of F{> 100 MeV) = 
1.4 X 10^^ phcm~^ s~^, corresponding to an average 
luminostity of (L-y) = 3 x lO^^ergs"^, and an aver- 
age photon index of (F) = 2.78 ± 0.03. In September 
2009, during a very bright flare reaching a peak flux 
of F(> 100 MeV) = 5.6 x 10"^ phcm'^ s'^, the hint 
of a spectral hardening accompanied the brightening 
of the source, possibly indicating a shift of the inverse 
Compton peak from 1-10 MeV (as typically observed 
in 3C 273) toward higher energies. 

When looking at the multi-band behaviour, the am- 
plitude of the variations seems to constantly increase 
from the X-rays up to the gamma-rays, with the MeV- 
GeV range showing the largest variations along the 
electromagnetic spectrum. On one side the correla- 
tion of the gamma-rays with the millimeter jet and 
with the optical flaring component [H, [l^, and on 
the other side the lack of correlation between the hard 
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X-rays and the millimeter emission [6| might indicate 
a different origin of the high-energy emission below 
and above the peak of the inverse Compton hump. 

An extensive correlation study and an accurate 
modelling of the broad band spectral energy distribu- 
tion during several epochs 20, 22] will allow to shed 
more light on the origin of the gamma-rays in 3C 273 
and their connection to the emission in other bands. 
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